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Abstract 

The cross section for e + e~ — > it + it~ J/tp between 3.8 GeV and 5.5 GeV is measured with a 967 fb _1 
data sample collected by the Belle detector at or near the T(nS) (n = 1, 2, 5) resonances. The 
y(4260) state is observed, and its resonance parameters are determined. In addition, an excess of it + it~ J/ip 
production around 4 GeV is observed. This feature can be described by a Breit-Wigner parameterization 
with properties that are consistent with the Y (4008) state that was previously reported by Belle. In a study 
of F(4260) — > -K + 7r~J/t/j decays, a structure is observed in the Mfa^J/ip) mass spectrum with 5.2a 
significance, with mass M = (3894.5 ± 6.6 ± 4.5) MeV/c 2 and width V = (63 ± 24 ± 26) MeV/c 2 , where 
the errors are statistical and systematic, respectively. This structure can be interpreted as a new charged 
charmonium-like state. 

PACS numbers: 14.40.Rt, 14.40.Pq, 13.66.Bc, 13.25.Gv 
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The F(4260) state was first observed by the BaBar Collaboration in the initial-state-radiation 
(ISR) process e + e~ —t 7isr7t + 7t~ J/ip []lj] and then confirmed by the CLEO [|2|] and Belle ex- 
periments ®\ using the same technique. Subsequently, a charged Z(4430) ± charmonium-like 
state was reported in the Tr ± ip(2S) invariant mass spectrum of B — > K7r ± ip(2S) [4] and two Z ± 
states were observed in the ir^Xci invariant mass distribution of B — > Kxd^ OU- Motivated 
by the striking observations of charged charmonium-like [0, [H] and bottomonium-like states |@], 
we investigate the existence of similar states as intermediate resonances in F(4260) — > n + n~J/ip 
decays. 

After the initial observations of the F(4260) fl-Hl, CLEO collected 13.2 pb" 1 of e+e~ data 
at yfs = 4.26 GeV and investigated 16 possible F(4260) decay modes with charmonium or light 
hadrons in the final state [7]. An ISR analysis by the Belle experiment with 548 fb _1 of data 
collected at or near yfs = 10.58 GeV [8] showed a significant y(4260) signal as well as an excess 
of n + n~J/ip event production near 4 GeV that could be described by a broad Breit-Wigner (BW) 
parameterization — the so-called F(4008) . Recently, the BaBar Collaboration reported an updated 
ISR analysis with 454 fb _1 of data and a modified approach for the background description [SO; 
the F(4260) state was observed with improved significance, but the F(4008) structure was not 
confirmed. Instead, they attributed the structure below the Y (4260) to exponentially falling non- 
resonant n + ix~ J/ip production. 

In this Letter, we report cross section measurements for e + e~ — > ix + ix~ J/ip between 3.8 GeV 
and 5.5 GeV, and a search for structures in the 7i + ir~ J/ip, ir^J/ip, and n + ix~ systems. The results 
are based on the full Belle data sample with an integrated luminosity of 967 fb 1 collected at or 
near the T(nS) resonances (n — 1, 2, 5). The Belle detector operated at the KEKB asymmetric- 
energy e + e collider [fioh and is described in detail elsewhere [Hill . We use the PHOKHARA 1 121 
program to generate signal Monte Carlo (MC) events and determine experimental efficiencies. 
The results reported here supersede those of Ref. 08J, wherein a subset of the Belle data sample 
was used. 

The event selection is described in Ref. B8fl. We require four well reconstructed charged tracks 
with zero net charge. For each charged track, a likelihood Cx is formed from different detector 
subsystems for particle hypothesis X e {e, fj,, it, K, p}. Tracks with a likelihood ratio TZk = 
c * c < 0.4 are identified as pions with an efficiency of about 95%. Similar ratios are also 
defined for lepton-pion discrimination [|13[|. For electrons from J/ip — >■ e + e~, one track should 
have 1Z e > 0.95 and the other track 1Z e > 0.05. For muons from J/ip — >■ /x + ^ , at least one 
track should have 71 ^ > 0.95; in cases where the other track has no muon identification, in order 
to suppress misidentified muon tracks, the polar angles of the two muon tracks in the fi + fi~ 
center-of-mass (CM) frame must satisfy | cos < 0.7. Events with 7 conversions are removed 
by requiring 7Z e < 0.75 for the n + ix~ candidate tracks. Furthermore, in J/ip — >■ e + e~, such 
events are further reduced by requiring the invariant mass of the ix + n~ candidate pair to be larger 
than 0.32 GeV/c 2 . Events with a total energy deposit in the electromagnetic calorimeter (ECL) 
above 9 GeV are removed in the J/ip — »■ e + e" mode because the MC simulation of the trigger 
efficiency for these Bhabha-like events does not accurately reproduce the data. There is only one 
combination of n + n~£ + £~ (£ = e, //) in each event after the above selections. 

Candidate ISR events are identified by the requirement |M r 2 ec | < 2.0 (GeV/c 2 ) 2 , where M 2 ec = 
(Pcm — P-K+ — P-k- — Pe+ — Pi-) 2 an d Pi represents the four- momentum of the corresponding parti- 
cle or composite in the e + e~ CM frame. Clear J/ip signals are observed in both the J/ip — >■ e + e~ 
and //+ it modes. We define the J/ip signal region as 3.06 GeV/c 2 < M{l + l~) < 3.14 GeV/c 2 
(the mass resolution for lepton pairs being about 20 MeV/c 2 ), and J/ip mass sidebands as 
2.91 GeV/c 2 < M(£+l~) < 3.03 GeV/c 2 or 3.17 GeV/c 2 < M (£+£-) < 3.29 GeV/c 2 , which are 
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three times as wide as the signal region. 

Figure (Ha) shows the 7T + 7r~£ + £~ invariant mass [14] distributions after all of these selection 
requirements are applied. Also shown in this figure are the background estimates evaluated using 
the normalized J/tp-mass sidebands. Two enhancements — the F(4008) and the F(4260) — 
above 3.8 GeV/c 2 are observed, consistent with the results of Ref. H but in disagreement with 
those of Ref. [9]. Other possible background sources not included in the sidebands are found to 
be small from MC simulation QV|] ; these include (1) n + n~ J /ip with J/ip decays into final states 
other than lepton pairs and (2) XJ/ip, with X not being a n + n~ pair, such as K + K~ or 7r + 7r~7r°. 
Non-ISR production of e + e~ — > TT + n~J/ip final states, such as e + e~ — > 77*7* — > 7p° J/ip, is also 
estimated to be small lfl5ll . FigureQlb) shows the measured cross sections for e + e~ — > 7r + 7r~ J/ip, 
where the error bars are statistical only. 
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FIG. 1: (a) Invariant mass distributions of it + tt~£ + £~. Points with error bars are data, and the shaded 
histograms are the normalized J ftp mass sidebands. The solid curves show the total best fit with two 
coherent resonances and contribution from background. The dashed curves are for solution I, while the dot- 
dashed curves are for solution EL The inset shows the distributions on a logarithmic vertical scale. The large 
peak around 3.686 GeV/c 2 is the ip(2S) — > 7r + 7r J/tp signal, (b) Cross section of e + e — > ir + ir J /if) 
after background subtraction. The errors are statistical only. 



Systematic uncertainties of the cross section measurement are found to be 7.9% and 7.3% for 
the e + e~ and yU + /i~ modes, respectively. The particle identification (PID) uncertainties, measured 
from pure 4>(2S) events in the same data sample, are 4.7% and 3.6% for the e + e~ and 
modes, respectively. Tracking efficiency uncertainties are estimated to be 3.3% for both e + e~ and 
modes in the momentum and angular regions of interest for signal events. The uncertainties 
associated with the choice of the J/tp mass window and |M r 2 ec | requirements are also estimated 
using pure ip(2S) events. It is found that MC efficiencies are higher than data by (4.5 ± 0.4)% in 
the e + e~ mode and (4.1 ± 0.2)% in the fi + mode. The differences in efficiencies are corrected 
and the uncertainties in the correction factors are incorporated into the systematic errors. Overall, 
together with the |M r 2 ec | requirements, these uncertainties contribute 0.6% for the e + e~ mode and 
0.3% for the mode within the J / ip mass window. Belle measures luminosity with a precision 
of 1.4% using wide-angle Bhabha events. The phokhara generator calculates ISR with 0.1% 
accuracy I112ll . The dominant uncertainties due to the MC generator are from three-body decay 
dynamics. MC simulation with modified invariant mass distributions weighted according to 
data distributions yields a 2% to 5% efficiency difference compared with a phase space n + ix~ mass 
spectrum model. Thus, we conservatively use 5% as the systematic error due to the approximations 
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made in the MC event generator. According to MC simulation, the offline trigger efficiency for 
four-track events of the studied topology exceeds 99%. A 1 .0% systematic error is included for the 
trigger uncertainty. The uncertainty of B(J/ip — > £ + £~) = B(J/ijj — > e + e~) + B(J/ip — > 
is 1.0% [fl^l. 

As a validation of our analysis, we also measure the ISR ift(2S) production rate using the same 
selection criteria. The cross sections are (14.12 ± 0.18 ± 0.85) pb and (15.13 ± 0.11 ± 0.79) pb 
at y/s = 10.58 GeV QJ/ZD for the e + e~ and modes, respectively. Our measurement agrees 

within errors with the prediction of (14.25 ± 0.26) pb II 1811 using the world-average resonance 
parameters 111 611 . 

An unbinned maximum likelihood fit is performed to the n + n~J/i/) mass spectrum above 
3.8 GeV/c 2 . As there are two enhancements observed, as shown in Fig. [TJaX we use the same 
fit strategy as in Ref. [81]. Two coherent BW functions (R 1 , R 2 ) are used to describe the F(4008) 
and F(4260) structures, assuming there is no continuum production of e + e~ — > n + ix~ J /ip. In the 
fit, the background term is fixed at the level obtained from a linear fit to the sideband data. The 
solid curves in Fig. [TJa) show the fit results. There are two solutions of equal optimum fit quality. 
The masses and widths of the resonances are the same for the two solutions; the partial widths to 
e + e~ and the relative phase between the two resonances are different (see Table U) [19]. The fit 
quality is estimated using the reduced \ 2 statistic; we obtain x 2 /ndf = 101.6/84, corresponding 
to a confidence level of 9.3%. Systematic uncertainties in the extracted values of the resonance 
parameters arise from the absolute energy scale, the detection efficiency, background estimation 
and the parameterization of the resonance models. The absolute energy uncertainty is estimated 
from the ip(2S) mass fit. Uncertainty in the detection efficiency does not affect the mass and width 
measurements, but could affect the measurement of the partial width to e + e~. Systematic uncer- 
tainties associated with the background contribution are estimated by varying the background level 
by ±lcr in the fit. Resonance parameterization is studied by changing the F(4260) BW function 
from a parameterization with a constant width to another with a three-body phase- space-dependent 
function. The interference between the two resonances, F(4260) and F(4008), depends on the 
structure of the n + n~J/?p amplitude, which can be different for the two resonances. We conserva- 
tively estimate possible systematic effects by performing a fit without the interference between the 
F(4008) and the F(4260) and taking the difference compared to results with interference as the 
corresponding systematic error. All of these contributions are summarized in Table U The mea- 
sured mass, width and the product T ee B(Ri — > n + n~J/iJj) are consistent within statistics with the 
previous results M. 12011 . 

Figure |2] shows the Dalitz plot for events in the F(4260) signal region (4.15 GeV/c 2 < 
M(n + it~ J/ip) < 4.45 GeV/c 2 ), where we observe structures in the n + n~ and n + J/ r tp systems. 
The inset is for the events in the J/^-mass sidebands, where no obvious structures are observed 
in the non-7r + 7r~ J /if) background events. 

Figure [3] shows a projection of the M(n + ix~), M(n + J/ip) and M(7i~J/ip) invariant mass 
distributions for events in the y(4260) signal region. Background contributions are estimated 
from the normalized J/ip mass sidebands. There are / (980), / (500) and non-resonant S-wave 
amplitudes in the n + n~ mass spectrum. In the n^J/ip mass spectrum, there is a significant peak 
around 3.9 GeV/c 2 (called the Z^OO)^ 1 hereafter) that we interpret as evidence for an exotic 
charmonium-like state decays into n ± J/?p; the broader peak near 3.5 GeV/c 2 is a reflection of the 
Z^OO^, as confirmed by MC simulation. The pure S-wave amplitudes that describe the 

7r + 7r~ invariant mass distribution well cannot reproduce the structure at 3.9 GeV/c 2 in the n ± J/ip 
mass spectra, as shown in Fig. |3] with the open histograms. 

An unbinned maximum likelihood fit is performed to the distribution of M max (irJ/ijj), the 
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TABLE I: Results of the fits to the 7r + 7r~ J/tp mass spectrum with two coherent resonances. M(Ri), 
r to t(^i) and T ee B(Ri ->■ it + it~ J/tp), i = 1, 2 represent the mass (in MeV/c 2 ), total width (in MeV/c 2 ) 
and product of the branching ratio for the decay into ir + ir~ J/ip and the e + e~ partial width (in eV/c 2 ) 
for the two resonances, respectively. The parameter 4> (in degrees) is the relative phase between the two 
resonances. The first and second errors are statistical and systematic, respectively. 



Parameters Solution I Solution II 

M{Ri) 3890.8 ±40.5 ± 11.5 

r tot (i?i) 254.5 ±39.5 ± 13.6 

T ee B{Ri -> ir + ir~ J/ij>) (3.8 ± 0.6 ± 0.4) (8.4 ± 1.2 ± 1.1) 

M(R 2 ) 4258.6 ± 8.3 ± 12.1 

r tot (i? 2 ) 134.1 ± 16.4 ± 5.5 

T ee B(R 2 -> vr+vr- J/ip) (6.4 ± 0.8 ± 0.6) (20.5 ± 1.4 ± 2.0) 

4> 59±17±11 -116 ± 6 ± 11 



2 



- 1.5 
1 






0.5 






: ..* : (.r. V-.;*.. s 

. *• -{«-•. s. s *••.;*. 
• .. -. V. "- • 






12 14 16 18 


• ' 

v »•.■* -'5''. ''* ■***'• 











3 10 12 14 16 18 20 22 



MVJ/y) (GeV/c 2 ) 2 



FIG. 2: Dalitz plot for y(4260) -»• vr+vr" J/^ decays for 4.15 GeV/c 2 < M(vr+7r- J/V>) < 4.45 GeV/c 2 . 
The inset shows background events from the J/ip-mass sidebands (not normalized). 




FIG. 3: Invariant mass distributions of (a) 7r + 7r~, (b) ir + J/ip and (c) ir~ J/ip for events in the Y(4260) 
signal region. Points with error bars represent data, shaded histograms are normalized background esti- 
mates from the J/^-mass sidebands, solid histograms represent MC simulations of 7r + 7r~ amplitudes B21I1 
(normalized J/^-mass sideband events added) and dashed histograms are MC simulation results for a 
Z(3900) ± signal 
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maximum of M(n + J/ip) and M(ir~ J/i[)). The signal shape is parameterized as an S-wave BW 
function convolved with a Gaussian whose mass resolution is fixed at the MC-estimated value of 
7.4 MeV, and the background is approximated by a cubic polynomial. The mass-dependent effi- 
ciency is also included in the fit. Figure |4]shows the fit results. The fit yields 159 ±49 ±7 Z(3900) ± 
events, with a mass of (3894.5 ± 6.6 ± 4.5) MeV/c 2 and a width of (63 ± 24 ± 26) MeV/c 2 , where 
the errors are statistical and systematic, respectively. The largest contributions to the systematic 
uncertainties arise from the parameterization of the signal and background shapes, the mass resolu- 
tion and mass calibrations. The statistical significance of the Z^OO^ state is found to be 5.5<r in 
the nominal fit, and is larger than 5.2a in all alternate fits used for systematic checks with different 
background shapes, fit ranges and BW resonant models. The significance is calculated by compar- 
ing the logarithmic likelihoods with and without the Z(3900) ± signal, and including the change of 
the number of parameters in the fits. From the signal yields and the MC- simulated efficiencies, we 
obtain the ratio of the production rates m^t^^^^m^^I^M = (29 .o ± 8 .9)%, 
where the error is statistical only [|22~ll . We test the hypothesis that interference between the S- 
and D-waves in the ix + n~ system might produce a structure similar to the enhancement observed 
in data. Although the statistics do not allow us to fully explore the Dalitz plot via an amplitude 
analysis, we find that these partial waves alone cannot produce an ± J/?p invariant mass peak near 
3.9 GeV/c 2 112 111 . Inclusion of a Z^OO)^ 1 with the mass and width determined above significantly 
improves the agreement between predicted and observed Dalitz plot distributions. 




FIG. 4: Unbinned maximum likelihood fit to the distribution of the M max (irJ/ip). Points with error bars 
are data, the curves are the best fit, the dashed histogram is the phase space distribution and the shaded 
histogram is the non-7r + 7r~ J/ip background estimated from the normalized J/ip sidebands. 

In summary, the cross section of e + e~ — > n + ix~ J /ip is measured from 3.8 GeV to 5.5 GeV. 
The F(4260) resonance is observed and its resonant parameters are determined. In addition, the 
F(4008) state is confirmed. The intermediate states in F(4260) — > n + ix~ J /ip decays are also 
investigated. A 2T(3900) ± state with a mass of (3894.5 ± 6.6 ± 4.5) MeV/c 2 and a width of 
(63±24±26) MeV/c 2 is observed in the ^J/ip mass spectrum with a statistical significance larger 
than 5.2a. This state is close to the DD* mass threshold; however, no enhancement is observed 
near the D*D* mass threshold. As the zT(3900) ± state has a strong coupling to charmonium and 
is charged, we conclude it cannot be a conventional cc state. 
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for efficient solenoid operations; and the KEK computer group, the Nil, and PNNL/EMSL for 
valuable computing and SINET4 network support. We acknowledge support from MEXT, JSPS 
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Note added- As we were preparing to submit this paper, we became aware of a paper from the 
BESIII Collaboration Q that also reports on the Z(3900) ± . 
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